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1.  Introduction 

When  atmospheric  air  Is  subjected  to  either  a  heating  source  or  an  external 
Ionizing  source,  the  chemical,  electrical,  and  thermodynamic  properties  of  the 
air  become  altered.  Of  particular  Interest  are  the  properties  of  air  plasmas 
formed  by  sources  such  as  charged  particle  beams,  beams  of  electranagnetlc  radia¬ 
tion,  and  electrical  discharges  passing  through  the  air.  Such  sources  can  Ini¬ 
tiate  a  complicated  and  diverse  set  of  reactions  which  radically  change  the  sub¬ 
sequent  composition  and  behavior  of  the  air.  A  computer  code  which  models  In 
detail  the  chemical  formation  and  decay  of  air  plasmas  would  thus  be  useful  for 
a  wide  range  of  applications. 

This  report  describes  an  air  chemistry  code,  CHHAIR,  developed  to  predict 
the  behavior  of  air  plasmas  generated  by  high  power  lasers.  The  code  follows  the 
time  histories  of  many  particle  species  and  several  temperatures  under  the  Influ¬ 
ence  of  external  heatinq/lonizing  sources.  CHMAIR  additionally  calculates  the 
radiation  emitted  over  a  wide  band  from  the  disturbed  air.  The  radiation  aspect 
of  hot  air  will  be  the  subject  of  another  report.  This  report  also  presents  rel¬ 
evant  electron  Impact,  excitation  and  Ionization  rate  coefficlerts,  based  on  a 
Maxwellian  electron  velocity  distribution,  for  the  species  of  interest. 

The  air  chemistry  code  contains  no  hydrodynemics;  however.  It  Is  contempla¬ 
ted  that  this  code  or  a  simpler  version  of  It  will  be  coupled  to  a  hydrodynamics 
code.  A  self  consistent  hydro  cooe  with  an  air  chemistry  should  be  of  conslder- 
able  Interest  In  the  understanding  of  air  plasmas. 

Note:  Manuacript  mbmitted  Aufuit  23.  1979. 
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2.  Disturbed  Air  Species  and  Reactions 


Pure  filr  at  sea  levels  for  our  purposes,  consists  mainly  of  and  0^.  When 
this  mlvture  Is  subj«cted  to  the  Influence  of  an  loni zing/heating  source,  posi¬ 
tive  Ions,  free  electrons,  negative  Ions,  atomic  species,  and  molecuUr  species 
consisting  of  various  coiii)1  nations  of  atomic  nitrogen  and  oxygen  are  formed. 

Thus,  our  first  generation  air  chemistry  code,  CHMAIR,  considers  the  time  his¬ 
tories  of  the  following  species:  ^2*  ^°2’  n'*’,  0^, 

0~,  Oj’  and  e.  In  addition  to  those  species,  the  code  calculates  the  time  his¬ 
tories  of  the  electron  temperature,  T^,  the  vibrational  temperature.  T^,  end 
the  kinetic  temperature,  T^,  of  the  heavy  particles. 

When  a  volume  of  air  Is  Ionized,  the  formation  and  the  decay  of  the  plasma 
are  controlled  by  a  large  nunt>er  of  atomic,  molecular  and  chemical  processes. 
These  processes  are  Indicated  below  In  the  details  of  the  reaction  kinetics  for 
the  specleii  of  Interest  to  the  current  air  chemistry  code,  CHnniR.  Ko««ver.  It 
should  be  stated  at  the  outset  that  several  minor  atmospheric  species,  especially 
CO^  and  H^O  could  affect  air  breakdown  and  also' alter  the  ensuing  air  chemistry. 
These  species  will  be  considered  In  the  development  of  the  second  ge>«rat1on 
CHNAIR  code^ 

Processes  Included  Reaction  Rate  (Coefficient  (cm  /sec)  Referef'ces 

Electron  Impact  lon.zatlon: 


e+N^^e  +  e  +  N* 

See  Table  1 

1 

e  +  K^e  +  e-frlf'" 

See  Table  1 

1 

e  +  O^^s  +  e  +  O* 

See  Table  1 

2 

2 


P^tssis  IncTtidtd 


e+O^e+e+0 


Electron  Iii|>act  Dissociative 
Ionization: 


e  +  H^-H+N+e 


e  +  Og-^O  +  O  +  e 


Radiative  Recombination: 


e  +  N  ♦  N  +  hi 


e  +  0  ♦  0  +  hi 


Three>Body  Recombination: 


e  +  e  +  N  ■►e  +  N 


e+e+0  *6+0 


e+A  +m*A+M 


Dissociative  Recombination: 


^2+6  -Jl  +  N 


O2  +  e  ♦  0  +  0 


NO  e  N  ^  0 


Reaction  Rate  (Coefficient  (cm'^/sec)  References 


See  Table  1 


See  Fig.  1 


See  Fig.  1 


2xlO’^\"°*^[  1.75- log 
75-log  T 

c  6 


.  -27  -a.«; 

8X10  ~  T  - 
e 


8x10-2'  ^^-4.5 
4xi0-’l  T/2-' 


.  -  ,--8  T  -0-39 
4.3x10  Tg 


1.5x10"®  T  <  0.1 

e  ’  e  “ 


2.1x10"®  T  T  >  0.1 
e  e 

2x10"® 


9.8xlO"®T 


Charge  Exchange: 

nJ  +  Og  ♦  N2  *  °2 

N2+  N  ♦  N2  +  N* 
N%  NO  ♦  N,  +  NO^ 

nJ  +  0  ♦  N2  +  o' 

O2  +  NO  -  O2  +  NO”^ 
+  0  >  N  +  O"^ 

N*  +  O2  ♦  N  +  O2 

+  NO  ♦  N  +  no'*’ 

o'*"  +  O2  ♦  0  +  0* 

O'*'  +  NO  ♦  0  +  NO"*^ 


RMCtloft  Rate  CCoefflcient  (cw^/tec)  Referencet 


2,7x10”*^  Tg  <  0.3 

4.2x10“'^  T„  >  0.3 

9  g 

10-12 

3.3x10-1° 

4.2x10-12  Tg-°*2^,  Tg  <  0.13 
1.6x10-11  T  °*^1,  T  >  0.13 

g  g 

4.4x10-1° 


10 


-12 


2.8x10-1°,  Tg  <  0.39 

4.8x10-1°  T  °*^1^.  T  >  0.39 
9  9 


8xlC 


-10 


4.6x10-12  T/°’^,  T„  <  0.155 
9  9  - 

lQ-10  1.2  >0.155 

S  S 

7.5  X  lO"!^,  Tg  <  0.1 

3.2x10-11  T  I’^'l,  T„  >  0.1 
9  9 


11 


12 

12,13 

14 


15 


e 

A 


2 


11  19 

A  A 


12 

11,16 


17,18 


4 


Lflllfp  i  Viu, 


3 

Pfoctst^s  IhC^jdtd  Reaction  Rate  (CoefficUnt  (tM  /stc)  Riftrinc^s 

Ion-Molecule  Rearrangements 

0"^  +  Ng  *  NO*  +  N  3.1  X  10’'^  Tg'^‘^.  Tg  <  0.065 

1.2x10"^°  T  0.065  <  <  0.67 

g  '  g  - 

8.5x10“^'  T  0.065  <  T.  <  0.67  11.16 

9  9 

N*  +  0  ♦  MO'^  +  N  2.Gxl0"^^  "  4.2xl0‘^^  Tg  <  0.13  H 

1.1x10“^°  Tg°*^  -  1.6x10“^^  ^ 


N*  +  O2  NO*  +  0 

2.8x10“^°,  Tg  <  0.39 

11 

ft 

4.8x10'^°  T  >  0.39 

9  ’  g 

• 

0*  +  N  ♦  no'*'  +  0 

.  .  ,„'10 
i.ZXlU 

IQ 
*  ^ 

Mutual  Neutralization: 

o'  ♦  H*  -  Oj  +  Nj 

2.5x10’®  Tg“®’^ 

7 

02"+  oj  ♦  02  +  Og 

f  6x10’®  T 

g 

7 

O2  +  NO*  ♦  O2  +  NO 

9.5x10’®  T  ’°*® 

g 

7 

0"  +  NO*  ♦  0  +  NO 

7.7x10’®  T  ’®*® 

g 

7 

p 

0  "+  N  *-  0  +  N 

4.1x10’®  T  ’®*® 

g 

7 

% 

0"  +  0*  ♦  0  +  O2 

1.6x10’®  T  ’°*® 

g 

7 

Proctists  loci  wild 

0“  +  0"^  ♦  0  +  0 

0~  +  NO*  ♦  0  V  NO 
Dissociative  Attactjment : 


Reaction  Rate  (Coefficient 


4.2xl0-»  T  -“-S 


7.7X10"*  T  ■*'* 


e  +  0  ♦  0"  +  0 
2 

Attachment: 

e  +  0  +  0  ♦O’  +  O 
2  2  2  2 

e  +  O  +  N  ♦O'  +  N 
2  2  2  2 

e  +  C  *  0"  +  hv 
Associctive  Detachment: 

0"  •»  0  ♦  0  +  e 
2 

0"  +  N  ♦  NO  +  e 

Colllslonal  Petachaent; 

O’  +  N  -0  +  N  +  e 
2  2  2  2 


See  Fig.  2 


3.5x10-31  yl  e-0‘^52/T^ 


10-31 


1.3x10 


-’5 


2x10' 


•10 


2x10 


-10 


4.8x10-^*^  T  e-°*^3/T 
9 


References 

7 

7 

20 

21 

21 

21 

21 

21 

21 
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fhrocMfis  Incl-ttdtd 

r.MCtlon  Rate  (Coefficient  (ce^/«fc) 

References 

Og  +  0^  ♦  O2  +  Og  e 

1.7x10"^  T 

9  ^ 

21 

Og  +  O^Og  +  O  +  e 

3.6x10"^^ 

21 

0"  +  Ng  ♦  0  Ng  +  e 

2.3x10"^ 

21 

0"  +  Og  •>  0  ■♦■  Og  +  e 

2.3x10"® 

21 

Therma]  Dissociation: 

+  Ng  -  N  +  N  ♦  Ng 

7.4X10-®  T  ^-9.75/T 

9  ^ 

22 

N2  +  O2  •»  N  +  N  +  O2 

2.2X10-®  T  -0-82 

9 

22 

N2  +  0-  N+  N  +  0 

2.2x10"®  T  "®*®2 

9  ^ 

22 

N2  +  N'^N*N-»-N 

-O  -1  K  _Q  7K/T 

4x10  «  T  *-  e  g 

9  ^ 

22 

O2  +  N2  ♦  0  +  0  +  N2 

2.7x10-10  T  -1-^  c-5-12/% 

9  ^ 

22 

02  +  02  -  0  +  0  +  02 

5.8x10"®  T  e"^‘12/l’g 

9  ^ 

22 

02  +  0-  0  +  0  +  0 

1.3x10"®  T  "1*0 

9 

22 

O2+NO-O+O+NO 

1.16xl0"10  T  "1*®  a"^‘l^^lg 

9  ^ 

22 

PFgmil.JigrlM^ 

RMctlon  Rate  (Coefficient  (cw^/iy;) 

References 

Og  +  N-  O  +  O  +  H 

1.16x10'^“  T 

9  * 

22 

NOg  +  H-  NO  +  O  +  H 

1.8x10“®  e"^*®*^^g 

22 

NO  +  N^^N+O  +  Nj 

1.5x10“®  e“®*®^^g 

22 

N0+02*N+0+02 

1.5x10“®  e‘®*®/^g 

22 

« 

N0♦N■►^«^0♦N 

1.05x10“®  T“^*®  e“®*®^^q 

9 

22 

N0  +  0*N  +  0  +  0 

1.05x10"®  T  e“®*®^^g 

9  * 

22 

NO-»-NO*N+0+NO 

1.05x10“®  T  e“®*^^^g 

9 

22 

Neutral  Rearrangeoent : 

N  +  Og  ♦  NO  +  0 

,  .^-10  ,  .-0.27/T 

i.j^xiu  1  e  •  q 

9  ^ 

12 

N(^D)  +  Og  ♦  NO  +  0 

4.7x10“^^  T 

9 

23,24 

0  +  Ng  -  NO  +  N 

1.3X10'!®  e“®'27/Tg 

12 

N  +  NO  ♦  Ng  +  0 

8.2x10"!! 

12 

N{^D)  +  MO  -  Ng  +  0 

6X10“!! 

12 

Processes  Included 

3 

Reaction  Rate  iCoefficlent  (cm  /sec) 

References 

0  -t  NO  N  +  O2 

-11  -1*6/T 

2.9x10  T  e  3 

g 

12 

N  +  NOg  *  NO  +  NO 

6x10"^^ 

25 

N  +  NO2  ♦  N2  +  O2 

6x10"’^ 

25 

N+N02-N2+0+0 

6x10"^^ 

25 

0  +  NO2  -  NO  +  O2 

i.7xl0‘^°  T 

9 

26 

NO  +  NO  +  O2  -  NO2  +  NOp 

,  0  in-39  ,0.046/T 

3.3x10  R  g 

26 

Neutral  Recombination: 

N  +  N  +  N2  -  N2  +  N2 

6.9x10'^^  T 

9 

22 

N+N+N-Np+N 

3.7x10“^^  T 

9 

22 

0  +  0  +  N2  O2  +  N2 

2x10“^^  T 

9 

22 

0  +  0  +  O2  +  O2  +  O2 

8.6x10"^^  T 

9 

22 

0  +  0  +  0*02  +  0 

1.9x10"^^  T 

9 

22 

O+O+N^Og+N 

1.7x10”^^  T 

9 

22 

9 


Processes  Included 


.3 


O  +  O  +  NO-Og  +  NO 

1.7x10'"  T 

9 

22 

N+0+N2-N0  +  N2 

1.1x10'^^ 

22 

N  +  0  +  02-N0  +  02 

1.1x10'^^ 

22 

N  +  N  +  O2  ♦  N2  +  Og 

2x10'^^  T 

9 

22 

N  +  N+  N0-H2  +  N0 

2x10'^^  T 

9 

22 

N+O  +  NO^NO  +  NO 

7.2x10'^^  T 

9 

22 

N  +  O  +  N-  NO+N 

7.2x10'^^  T 

9 

22 

N  +  O  +  O^NO  +  0 

7.2x10"^^  " 

9 

22 

0  +  NO  +  N2  ♦  NO2  +  N2 

1.5x10'^^ 

26 

0  +  NO  +  O2  NO2  +  O2 

3x10-33  ,^.08/Tg 

26 

Electron  Impact  Dissociation: 

e  +  0,  -  O  +  O  +  e  l,3xi0"^°  T  °’^(4.4+T  )e"®*^^^e  2,27 

+1.7xl0‘®  T  ^'^(2.25+T  )e"^’®^^e 
e  e 


10 


I  i  -If  lll'l'tfl 


Processes  Included  Reaction  Rate  (Coefficient  {on^/sec)  References 

e  +  Ng-N  +  N  +  e  6.2xlO"^°Tg°-®(4.8+Tg)e“^*^^^%  28 

Associative  Ionization; 

N  +  0  ♦  NO^  +  e  Detailed  balance  2 

N(^0)  +  0  ♦  NO*  +  e  Detailed  balance  2 

Photo-Detachment ; 

hv  +  O2  "•  O2  e  Wavelength  Dependent  29 

hv+0~-*0  +  e  Wavelength  Dependent  29 

♦Note:  Three  body  reaction  rata  coefficients  are  In  units  of  cm®/sec. 
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Table  I 


Electron  Impact  Ionization  Rate  Coefficients 


Te(aV) 

N2  loniz 

N  loftiz 

0  loniz 

0^  loniz 

0.2 

1.92  (-43) 

O.C 

3.14  (-38 

2.88  (-38) 

0.3 

6.14  (-32) 

0.0 

2.C8  (-28 

1.66  (-28) 

0.4 

3.65  (-26) 

0.0 

1.71  (-23: 

1.32  (-23) 

0.5 

1.09  (-22) 

7.24  (-23) 

i.55  (-20; 

1.18  (-20) 

0.6 

2.32  (-20) 

3.14  (-20) 

1.46  (-18] 

1.12  (-18) 

0.7 

1.08  (-18) 

1.61  (-18) 

3.81  (-17 

2.96  (-17) 

0.8 

1.96  -17 

2.88  (-17) 

4.43  (-16 

3.49  (-16 

0.9 

1.83  (-16) 

2.66  (-16) 

3.01  t-15) 

2.40  (-15) 

1.0 

1.15  (-15) 

1.57  (-15 

1.40  (-14; 

1.12  (-14 

1.1 

5.13  (-15) 

6.69  (-15) 

4.95  (-14 

4.05  (-14) 

1.3 

5.17  (-14) 

6.26  (-14) 

3.50  (-13] 

2.91  (-13) 

1.5 

2.86  (-13) 

3.25  (-13) 

1.49  (-12 

1.25  (-12) 

1.7 

1.07  (-12 

1.15  (-12) 

4.54  (-12 

3.85  (-12) 

1.9 

3.07  (-12) 

3.16  (-12) 

1.11  (-11] 

9.36  (-12) 

2.1 

7.25  (-12) 

7.2  (-12) 

2.29  (-11] 

1.94  (-11) 

2.3 

1.49  (-11) 

1.43  (-11) 

4.20  (-11] 

3.54  (-11) 

2.5 

2.74  (-11) 

2.55  (-11) 

7.03  (-11] 

5.83  (-11) 

2.7 

4.62  (-11) 

4.20  (-11) 

1.09  (-10] 

9.00  (-11) 

3.0 

8.99  -11) 

7.89  (-11) 

1.92  (-10] 

1.55  (-10) 

3.3 

1.56  (-11) 

1.33  (-10) 

3.05  (-10] 

2.44  (-10) 

3.5 

2.15  (-10) 

1.8  (-10) 

4.00  (-10 

3.14  (-10) 

3.8 

3.28  (-10) 

2.68  (-10) 

5.70  (-10 

4.41  (-10) 

4.1 

4.73  (-10) 

3.78  (-10) 

7.75  (-10] 

5.86  (-10) 

4.5 

7.21  (-10) 

5.60  (-10) 

1.10  (-9) 

8.10  (-10; 

5.0 

1.12  (-9) 

8.45  (-10) 

1.59  (-9) 

1.12  (-9) 

3.  Metastable  States 

In  Section  2  we  described  the  species  and  the  relevant  reactions  in  our 
multispecies  air  chemistry  code.  These  species  are  considered  to  be  in  their 
ground  states.  Howeve-,  several  metastable  states  arise  in  air  when  it  is  sub¬ 
jected  to  an  ionization  impulse.  These  are;  0(^0),  0(^S),  N(^D),  N^^P),  N2(A^i), 
02(a^&),  02(b^£),  O^(^O),  O^(^P),  N^(^O),  and  N^(^S).  Their  influence  on  emis¬ 
sion  and  deionization  processes  In  the  E  and  F  regions  of  the  ionosphere  have 

been  considered  in  detail.^*  30,  31  metastable  states  have  also 

31-33 

been  considered’  in  multispecies  codes  for  the  disturbed  D-region. 
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These  metastable  states  store  a  certain  amount  of  the  energy-  deposited  In 
the  disturbed  air.  Part  of  this  stored  energy  is  released  as  radiation  while  the 
major  part  is  converted  into  the  kinetic  energies  of  the  electron  gas  and  the 
heavy  particles.  This  is  due  to  deexcitation,  quenching  and  charge  exchange  pro¬ 
cesses  whose  rates  increase  with  increasing  electron  and  neutral  densities.  The 
emitted  radiation  from  these  and  short  lived  states  can  be  utilized,  in  principle, 
for  diagnostics  and  tracking  ot  the  disturbance.  In  addition  to  their  "influence 
on  the  kinetic  temperatures,  some  of  the  metastables  play  important  roles  in  the 
production  of  new  species.  Thus,  the  excitation  and  deexcitations  of  the  meta- 
stable  states  constitute  some  of  the  basic  elements  of  the  air  chemistry,  De- 
tailed  excitation  and  deexcitation  mechanisms  of  these  states  have  been  reviewed;’’^*'^^ 
in  this  section  we  present  the  important  processes  and  give  the  relevant  reaction 
rate  coefficients  utilized  in  the  air  chemistry  code  CHIV.IR: 

2 

(a)  NjfA  r):  This  state  is  one  of  the  lowest  triplet  states  of  N,  with  a 
threshold  excitation  energy  of  6.2  eV  and  a  lifetime"^  of  ~  1  sec.  It  is  ex¬ 
cited  from  the  ground  state  of  by  electron  impact  and  by  cascade  from  higher 

3  3 

electronic  states  of  N^,  especially  C  »  and  B  w,  which  are  also  excited  from  the 
ground  state  by  electron  impact.  The  electron  impact  excitation  rate  coeffi- 
cients  ’  of  the  triplet  states,  C  w,  B  w  and  A  z  are  given  in  Table  U  as  a  func¬ 
tion  of  the  electron  temperature.  It  should  be  stated,  that  the  excitation  of 
3  38 

C  T  State  is  one  of  the  strongest  for  air  or  pure  nitrogen  when  they  are  bom- 

3 

barded  by  high  energy  electron  beams.  The  A  i  state  is  deexcitated  by  electrons 
through  the  superelastic  collisions  and  quenched  predominantly  by  0^,  N,  and  0 
with  rate  coefficients  of  2.5x10"^^  cm^/sec  (Ref.  39),  5xl0"^^  cm^/sec  (Ref.  40) 
and  5x10  cm^/sec  (Ref.  41),  respectively. 
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Table  II 


Electron  Impact  Excitation  Rate  Coefficients 


C^i  a(^A) 


0.2 

1.3(-23) 

2.8< 

;-25] 

4.71 

;-33) 

5.1(-13' 

1.61 

[-14) 

0.3 

1.3(-18 

6.81 

-20 

8.0< 

-25) 

3.7(-12) 

3.2! 

!-13) 

0.4 

4.3(-16) 

3.51 

1.1( 

[-20) 

1.2(-11 

1.7! 

[-12) 

0.5 

1.5(-14] 

1.31 

1-15] 

3.21 

[-18) 

2.3(-ll) 

4.11 

[-12) 

0.6 

J.5(-13) 

1.9( 

:-14] 

1.4! 

1-16) 

3.8(-ll) 

7.91 

[-12) 

0.7 

8.3(-13 

1.2( 

;-i3] 

2.3( 

[-15 

5.5-11 

1.41 

[-11) 

0.9 

8.0(-12j 

1.4! 

1-12] 

7.91 

[-14) 

l.l(-lO) 

2.41 

-11) 

1.0 

1.8^11 

3.41 

:-i2] 

2.8| 

[-13) 

1.3-10) 

3.21 

(-11) 

1.2 

5.9(-llj 

1.31 

Ml] 

1.9| 

[-12) 

1.9(-10) 

4.61 

-11 

1.3 

9.4(-ll] 

2.1! 

[-11 

3.81 

[-12) 

2.4(-10) 

5.21 

(-11) 

1.5 

2.o(-io; 

4.9; 

[-11] 

1.21 

(-11) 

3.0{-10) 

6.71 

(-11) 

1.7 

?.4(-10] 

9.31 

[-11] 

2.91 

(-11) 

3.5(-10) 

7.91 

(-11) 

2.0 

4.3f-io; 

1  1.91 

[-10; 

7.8! 

(-11) 

4.4(-10) 

1.11 

(-10) 

2.5 

1.3(  -9) 

4.4! 

[-10] 

2.3! 

(-10) 

5.5(-10) 

1.31 

(-10) 

3.0 

2.1(-9) 

7.81 

[-10: 

4.81 

-10) 

6.5(-10 

1.51 

-10 

4.1 

3.7(-9) 

1.61 

t-9) 

1.21 

(-9) 

7.8(-10) 

1.81 

-10) 

5.0 

4.8(-9) 

2.41 

(-9) 

1.91 

(-9) 

8.5(-10) 

2.01 

(-10) 

3 

The  A  £  state  can  also  be  excited  through  the  resonant  charge  transfer 


M*  +  NO  ♦  Ng  (A)  +  NO'^ 


with  a  rate  coefficient  of  3x10“^®  cni^/sec. 

(b)  02(3^4)  and  02(b^i);  These  two  low  lying  metastable  states  of  O2 

36  42 

have  excitation  energies  of  ~  1  eV  and  1.63  eV,  respectively,  and  lifetimes 

of  '  4000  sec  and  12  sec,  respectively.  They  are  excited  predominantly  from  the 

ground  state  of  O2  by  electron  impact.  The  rate  coefficients  for  their  excita- 

t Ions  *  are  given  In  Table  II.  These  states  are  deexcitated  by  superelastic 

collisions  with  free  electrons.  Their  quenching  by  atmospheric  species,  however, 

1  1 

Is  not  fast.  The  effective  quencher  for  02(a  4)  Is  N  and  for  02(8  t)  1s  N2* 
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The  corresponding  rate  coefficients  are  3x10"'®  cm^/sec  (Ref.  #3)  and  2.5x10“ 
3  1 

cm  /sec  (Ref.  44),  respectively.  The  a  a  also  disappears  through  collisional 
detachment 

®2~  ♦  2O2  +  e 

and  associative  detachment 


0"  +  0,(3^)  ♦  0,  +  e 


with  rate  coefficients^®  of  2x10“'^  cm^/sec  and  3xl0"'®  cm^/sec,  respectively. 

2  2 

(c)  N(  D)  and  N(  P):  These  two  low  lying  metastable  states  of  N  have  exci¬ 
tation  energies^®  of  2.37  eV  and  3.56  eV,  respectively.  Their  lifetimes  are^® 

A 

6x10  sec  and  13  sec,  respectively:  They  are  excited  from  the  ground  state  of  N 

2  47 

by  electron  impact  and  corresponding  excitation  *  rate  coefficients  are  presen- 

2 

ted  in  Table  III.  N(  0)  is  excited  by  the  dissociative  recombinations 


nJ  +  e  -  N(^D)  +  N(^  S) 


and 


NO"^  +  e  -  N(^D)  +  0 

where  the  yield  of  N(^d)  in  the  first  reaction  is  near  unity^®  and  for  the  second 
40 

reaction  is  ~  0.76. 

The  low  lying  states  of  N  can  also  be  excited  through  electron  impact  disso¬ 
ciation  and  dissociative  ionization  of  N^.  However,  no  quantitative  measure¬ 
ments  for  the  yields  exist. 
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TiMte  1cN<  lyinf  states  are  deencltated  by  svperelastic  collisions  with  the 
free  electrons.  The  nost  leportant  lots  Mechanise  for  N<^)  i%  the  foreatlon  of 
NO,  via 

N(^D)  +  Og  NO  +  0. 

It  Is  also  quenched  by  the  atmospheric  specl-^s  N^,  0  and  NO  with  rate  coeffl* 
clents^^*’^  of  1.6x10“^*  cm^/sec.  1.8x10"^^  cm^/sec  and  7x10“^^  cm^/sec, 
respectively. 

Table  III 

Excitation  Rate  Coefficients  for  the  Low  Lying 
States  of  N 


T^(eV) 

^S-  2d 

S  -  p 

0.1 

8.0(-20) 

2.fi(*25) 

2.10(-14) 

0.2 

1.54(-14) 

2.03{-17) 

1.16(-11) 

0.3 

1.25(-12) 

1.11(-14) 

9.78(-ll) 

0.5 

4.2(-ll) 

1.47(-12) 

6.0(-10) 

0.7 

1.98(-10) 

1.63(-11) 

1.16(-9) 

1.0 

6.38(-10) 

8.91(-11) 

2.05(-9) 

1.2 

1.08{-9) 

1.73(-10) 

2.56(-9) 

1.5 

1.52(-9) 

3.38(-10) 

3.24(-9) 

2.0 

2.27(-9) 

6.54(-10) 

4.1(-9) 

3.0 

3.31{-9) 

1.23(-9) 

5.18(-9) 

5.0 

4.50(-9) 

1.99(-9) 

e.2(-9) 
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(d)  b{^D)  and  0(^S):  These  two  low  lying  metastable  states  of  0  have  exci¬ 
tation^  energies  of  i.96  eV  and  4.18  eV,  >*espect1vely  a.^d  lifetimes^  of  140  sec 

and  0  *ec,  respect Ively.  They  are  excited  from  the  ground  states  of  0  by  elec- 

2  47 

trc  1  ^ujpact.  The  corr“fespond1  ng  excitation  rate  coefficients  *  sre  given  In 
Table  IV.  They  are  also  excited  through  the  dissociative  recombination  of  0^. 

The  product  results^®  In  1.0  0(^P).  0.9  0(^0)  and  0.1  0('s).  respectively. 
Furthermore,  they  can  be  excited  by  electron  Inpact  dissociation.  However,  no 

accurate  cross  sections  for  those  processes  are  available.  Only  approximate  rate 

2  27 

coefficients  can  be  obtained  using  approximate  cross  sections. 

The  deexc^tation  of  these  states  proceeds  through  the  superelastic  colli¬ 
sions  with  the  free  electrons.  They  are  quenched  by  the  atmospheric  species  very 
effectively  at  atmospheric  densities. 


Table  IV 


Excitation  Rate  Coefficients  for  the  Low  Lying  States 

of  0 


T^(eV) 

3p  -  ^0 

3n  Is 

^0  -  h 

0.1 

1.92  (-18) 

1.78  (-28) 

2.25  (-19) 

0.2 

5.28  (-14) 

1.96  (-19) 

1.76  (-14) 

0.3 

1.76  (-12) 

2.10  (-16) 

8.07  (-13) 

0.5 

3.28  (-11) 

6.04  (-14) 

1.52  (-11) 

0.7 

1.21  (-10) 

7.25  (-13) 

5.4  (-11) 

1.0 

3.43  (-10) 

4.93  (-12) 

1.38  (-10) 

1.2 

5.20  (-10) 

1.06  (-11) 

1.97  (-10) 

1.5 

7.94  (-10) 

2.32  (-11) 

2.76  (-10) 

2.0 

1.21  (-9) 

5.15  (-11) 

3.98  (-10) 

3.0 

1.84  (-9) 

1.16  (-10) 

5.30  (-10) 

5.0 

2.52  (-9) 

2.21  (-10) 

7.16  (-10) 
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The  racst  cff«ctive  quencher$  of  0(^o)  are  N^.  0^  and  NO  and  their  q^enyM^i 
coefflclents^^*^®  Of  are  5x10“^^  cm^/sec.  7x10"^^  cm^/sec,  and  1.7x10“^^  cfi^/sec^ 
respectively.  For  0(^S),  the  most  effective  quencher  Is  0(^P)  w1?ere  the  rate  cp- 
efflclent^^’^^  Ts  8x10"^^  cm^/sec. 

(e)  0*(^D)  and  0*(^F):  These  low  lying  metastabie  states  of  0*'have  eititta- 
tlon  energies^®  of  3,32  eV  and  5.0  eV  and  lifetimes^  of  “  10^  sec  and  $  sllc,  if^^- 
spectlvely.  they  are  excited  from  the  ground  states  of  0*  by  low  energy^electrbrt 
Impact.  The  excitation  rate  coefficients  for  threshold  and  near  threshold  eitcl^' 
tatlon  derlved^^  from  the  relevant  collision  strengths  are  given  In  TabTe  V. 


These  states  also  arise  as  a  result  of  the  dissociative  Ionization  d^'O^. 

■  Cl  . .  •-  . 

However,  the  only  quantitative  data  on  this  aspect  Is  given  by  StebbIngS,'  fir' 
Incident  electron  energies  from  threshold  to  100  eV.  For  100  eV  el(Ktrons 

nn  O  ♦ha  nf  Oi^  efa>AC  1r5  r»<A  f)  and  0  l^Pl. 

W  •  »»  «•  ^  -m  p  v*..ar  -  %  V,  ..r..  ,  ,  ^ 

could  be  as  hlgh®^  as  0.44.  The  deexcitation  of  these  states  proceeds  through 

the  superelastic,  collisions  with  the  free  electrons.  However,  the  most  Important 

loss  mechanism  of  O^(^O)  Is  the  near  resonance  charge  exchange  with  Ng  and  0^ 

52  -9  3 

which  may  have  a  rate  coefficient  as  high  as  3x10  cm  /sec.  The  Importance 
of  these  reactions  Is  that  they  convert  an  atomic  Ion  Into  a  molecular  Ion  which 
has  a  larger  recombination  (dissociative)  rate  with  the  free  electrons. 


(f)  N'*’(^D)  and  n'‘^(^S):  These  low  lying  metastabie  states  have  excitation 
energies^®  of  1.89  eV  and  4.04  eV  and  lifetimes  of  250  sec  and  "■  1  sec,  resoec- 
tively.  They  are  excited  from  the  ground  state  of  N*  by  low  energy  electron  Im¬ 
pact.  The  rate  coefficients  for  these  excitations  based  on  threshold  collision 
strengths^^  are  given  In  fable  V.  These  metastabie  states.  In  principle,  should 
arise  through  the  electron  Impact  dissociative  Ionization  of  N  In  a  manner  slml 
lar  to  the  excitation  of  0^  metastabie  states  (through  the  dissociative  Ioniza¬ 
tion  of  0^  discussed  In  section  3e). 
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Transition 

0* 


Tahl*  V 

Exc1t«t<on  Rate  Coefficients  for  the  Low  Lying 
Netaetable  States  of  0^  and  N''^  Te  In  eV. 

Rate  Coefficient  (cw^/sec) 
3.2x10”®  (Tg)~^  ®  exp  (-3.33/T  ^ 

9.5x10“®  (r )”®*®  exp(-5.0/T.) 

6  a 

1.4x10”®  (r)”®*®  exp(-5.0AJ 

6  Q 

2.7x10”®  (r)”®*®  exp(-1.89/T  ) 
3.5x10"®  (r )"®*®  exp(-4.04/T  ) 

C  V 

6.5x10”®  (Tg)"®-®  exp(=2.15/Tg) 


The  deexcltatlon  of  these  states  proceeds  through  the  superelastic  collisions 
with  tnt  low  enerey  electrons.  Furthenaore,  the  charge  exchange  process  of  the 
Metastable  states  with  are  exothermic  and  should  convert  an  atonic  Ion  Into  a 
molecular  one.  However,  no  quantitative  data  exists  In  this  at'ea. 

The  table  referred  In  this  section  gives  the  excitation  rate  coefficients 
as  a  function  of  the  electron  temperature.  The  corresponding  deexcltatlon  rate 
coefficients  can  be  obtained  through  the  principle  of  the  detailed  balance. 

Finally,  the  first  generation  air  chemistry  code.  CHMRIR.  calculates  most  of 
these  netastable  states  through  steady-state  relations.  Their  effects  and 
their  reactions  are  thus  perpetuated  accordingly. 
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4.  Inelastic  Processes  and  the  Kinetic  Temperatures. 

Describing  the  development  and  decay  of  an  Inntzr^  gaseous  medium  r*equ1res  a 
good  knowledge  of  the  particle  Kinetic  energies.  This  Is  because  the  Ionization 
«ind  deionization  processes  depend  strongly  on  these  energies.  The  most,  accurate 
method  for  determining  these  energies  and  the  resultant  ionization  Is  to  solve 
the  appropriate  Boltzmann  c<)uat1ons.  coupled  to  the  particle  density  equations. 
Such  solutions  are  presently  Intractable,  however,  except  where  numerous  simplifi¬ 
cations  can  be  made.®^’^^  An  alternative  approach,  used  in  CHMAIR,  Is  to  use  Max¬ 
wellian  velocity  distributions  for  the  electrons  and  heavy  particles,  to  express 
the  collision  rates  as  functions  of  the  electron  and  heavy-particle  temperatures. 
The  validity  of  this  cssumption  has  been  extensively  studied  and  shown  to  be  not 
unreasonable  for  broakdoim  problems.  The  assumption  leads  to  an  equation  for 
the  electron  temperature  and  an  equation  for  the  heavy-particle  temperature. 

'’’he  processes  leading  to  energy  losses  by  electrons  In  are  Ionization, 

dissociative  ionization,  dissociation,  and  excitation  of  the  following  electronic 
3  3  3  3  1 

states;  A  i,  B  w,  C  w,  E  i,  a  w,  two  higher  electronic  groups,  and  the  vibra- 

12  3 

tional  excitations  of  the  ground  state.  The  relevant  rate  coefficients  ’’  for 
ionization,  dissociative  ionization,  and  excitation  of  the  electronic  states  are 
given  in  Tables  I-III,  VI.  The  electron  impact  excitation  rate  coefficients^  for 

55 

the  eight  ground  state  vibrational  levels,  based  on  measured  cross  sections^ 
are  given  In  Table  VII. 

For  Ng,  the  electronic  states  considered  as  energy  sinks  for  electrons,  are 

N^CB)  and  N^CA)  states.  The  excitation  rate  coefficients^  for  N2(B),  excited 

from  the  gn  md  state  of  based  on  a  measured  cross  sectlon^^,  are  gfven  In 

Table  VI.  As  for  the  excitation  of  N2(A)  no  measured  cross  section  Is  available; 

58 

however,  an  approximate  relation  based  on  Seaton  type  cross  section  with  a  gaunt 
factor  of  0.2  is  utilized. 
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T«b1«  VI 


Eltctroh  Ii^iact  Excitation  Rate  Coefficients 


E^t 

al. 

Opt. 

2  & 

Group 

4 

nJ(b) 

0.2 

2.0  (-35) 

7.0  (-31) 

1.6  (-14) 

0.3 

1.1  (-26) 

1.6  (-23) 

2.4  (-12) 

0.5 

1.0  (-20) 

1.3  (-17) 

3.2 

(-20) 

1.2  (-10) 

8.0 

(-21) 

0.7 

1.0  (-16) 

5  0  (-15) 

4.4 

(-17) 

7.2  (-10) 

1.1 

(-17) 

1.0 

1.7  (-14) 

2.5  (-13) 

1.0 

(-14) 

2.0  (-9) 

2.5 

(-15) 

1.2 

1.3  (-13) 

2.4  (-12) 

3.9 

(-13) 

3.2  (-9) 

1.0 

(-13) 

1.5 

9.0  (-13) 

1.4  (-11) 

8.4 

(-13) 

5.2  (-9) 

2.1 

(-13) 

2.0 

6.2  (-12) 

9.9  (-11) 

8.0 

(-12) 

8.3  (-9) 

2.0 

(-12) 

3.0 

4.0  (-11) 

4.6  (-10) 

8.4 

(-11) 

1.2  (-8) 

2.1 

(-11) 

5.0 

1.5  (-10) 

1.9  (-9) 

7.2 

(-10) 

1.3  (-8) 

1.8 

i-io) 
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ViDratlonal  Excitation  Rate  Coefficients 


CM 

X 


t  I 


I  I 


oToO  CJ^CM  ^tO^OCM^^rOOOCMCMCSiCVl  ^  ^ 

«t‘^00«~4mf-4tOOCMiMOOOO*^«Mf-l«M«-^rHOCMCsi^r<vf-4 

•  €•••••>  ••«•••••••••••%•• 


l-H  ^  I— <  f-H  1^  O' 

•  I  •  I  I  I  •  I 


o  « 2! 

T-l  4J  *S 
•  « 

'gr 

t- 

c 
M  at 


>0  o 


0'0vo0«no»«»<-ifn^»n»ot^r's.r^r«.r^r^i»si/>c»jevj<-«io»o 

^  ^  ^  f-^  eH  ^  ^  eM  0\  00 


<->4CMco^kn^r<<'COO>o<-iCMco^m«or^ooo^inoinOir>o 

•  ••••••••••eee«««<.«eeec»ee 

O  00000000«-*4  ^  ^  ^  ^  ^  ^  0-4  ^  f-'i  CM  to  ‘RT  ^  i/> 
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^atiui  -.'rTt^  .^-■.4- ,  an,‘  t] 


■i 

J 


The  Inelastic  pf*ocesses  through  which  electrons  lose  energ^^  In  N  arc  loniza- 

2  2 

tion  and  the  excitation  of  the  electronic  states  D,  P,  and  the  optically 

allowed  transitions,  especially  those  with  large  oscillator  strength.  The  rele- 

1  47 

vant  rate  coefficients  *  are  given  In  Tables  I  and  III.  As  for  N  ,  the  excita- 

J-  1  X  1 

tions  of  N(  p),  N  (  S)  and  the  optically  allowed  transitions  constitute  the  main 
Inelastic  processes  for  energy  loss  by  electrons.  The  relevant  ri.te  coefficients 
for  ^0  and  ar?  given  In  Table  V, 

The  Inelastic  energy  loss  processes  for  electrons  In  0^  are  Ionization,  dis¬ 
sociative  Ionization,  and  the  excitation  of  the  electronic  states,  a^a,  b^r,  A''z, 

3 

and  B  z  where  the  last  two  excitations  are  considered  to  lead  to  the  dissociation 

1  2 

of  0^.  The  relevant  rate  coefficients  *  are  given  In  Section  II.  As  for  the 
excitations  of  the  electronic  states  of  O^,  no  relevant  experimental  cross  sec¬ 
tions  are  available,  and  thus  one  can  only  provide  estimates  for  their  rate  coef- 
ficlents  based  on  Seaton  type  cross  sections. 

The  electron  energy  loss  processes  in  0  considered  in  CHMAIR  are  ionization 
and  the  excitation  of  the  electronic  states  ^D,  and  the  optically  allowed  tran¬ 
sitions.  The  relevant  rate  coefficients  are  given  in  Tables  1  and  IV.  As  for  0^, 
2  2 

the  excitation  of  0  and  P  and  the  optically  allowed  transitions  constitute  the 
main  processes  for  the  ewirgy  loss  by  electrons.  The  relevant  rate  coefficients 
are  given  In  Table  V, 

In  addition  to  the  above  Inelastic  processes,  the  electrons  lose  their 
energy  also  through  elastic  collisions  with  the  neutrals.  Coulomb  collisions  with 
the  Ions,  and  by  Bremsstrahlung  radiation.  Therefore,  the  time  history  of  the 
electron  temperature  Is  described  in  detail,  given  the  source  of  the  original 
heating,  the  detailed  loss  processes  discussed  In  this  section,  and  the  corre¬ 
sponding  deexcitation  processes  which  heat  the  low  energy  electrons. 


The  air  chemistry  code  calculates  two  other  temperatures,  the  vibrational 
temperature,  and  the  heavy  particle  ten^erature.  The  vibrational  temperature 
arises  as  a  result  of  electron  impact  collisions  with  the  ground  state  of  the  ni¬ 
trogen  molecule  where  eight  vibrational  levels  are  considered  to  be  In  a  vibra¬ 
tional  equilibrium.  The  electron  Impact  excitation  rate  coefficients  for  these 
levels  are  given  In  Table  VII,  and  the  corresponding  deexcitation  rates  are  obtained 
through  the  principle  of  detailed  balance.  Several  other  processes  also  result 
In  the  vibrational  excitation  of  Ng.  One  of  these  Is  the  quenching  of  0(^0)  by 
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N2  where  30X  of  the  Internal  energy  Is  converted  Into  the  vibrational  energy  of 
N^.  The  atom-atom  Interchange  forming  N2 

N  +  NO  ♦  N2  +  0 

also  contributes  to  the  vibrational  temperature,  where  one  quarter®^  of  the  exo¬ 
thermic  energy  Is  expended  Into  vibrational  energy.  The  radiative  cascade  from 
the  N2(A  i)  state  also  raises  the  vibrational  temperature.  Processes  which  result 
In  the  loss  of  vibrational  energy  are  the  deexcitation  of  the  vibrational  states 
by  electrons  through  the  superelastic  collision  and  quenching  by  neutral  species. 

The  heavy  particle  temperature  arises  as  a  result  of  dissociative  recombina¬ 
tion,  electron  collisions,  charge  exchange  and  neutral  rearrangement  and  quench¬ 
ing  processes  where  some  of  the  exothermic  energies  are  released  as  kinetic.  The 
dissociative  recombinations  are  an  nportant  source  of  neutral  kinetic  tempera¬ 
ture  where  the  potential  energy  Is  converted,  partially,  to  kinetic  energy. 

Other  processes  affecting  the  neutral  kinetics  temperature  are  the  vibrational- 
translational  Interchange. 
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5.  Reroarks  on  the  Reaction  Rate  Coefficients 

It  is  obvious  from  this  report  that  a  very  large  number  of  reactions  occur  in 
an  air  plasma.  Furthermore,  almost  all  of  these  reactions  depend  on  one  tempera¬ 
ture  or  another,  which  makes  the  kinetic  temperatures  of  the  plasma  important  ele¬ 
ments  of  the  deionization  processes  in  air.  However,  a  large  number  of  these 
reactions  have  been  measured  only  at  room  temperature,  or  slightly  higher.  There¬ 
fore,  the  extension  of  these  rates  to  higher  temperatures  may  be  in  arror.  Some 
reactions  have  been  measured  over  a  wide  temperature  range  and  yet  with  different 
temperature  dependencies.  A  case  in  point  is  the  dissociative  recombination  of 

NO^  where  one  measurement^^  gives  a  temperature  dependence  of  while 

9  -0  83 

another  measurement  gives  a  temperature  dependence  *  .  These  two  rates 

are  illustrated  in  Fig.  3  where  one  recombination  rate  coefficient  is  faster  chan 
the  other  by  a  factor  of  3  or  more  at  higher  temperatures.  Thus  the  inverse  pro¬ 
cess  of  associative  ionization  will  clearly  be  different,  depending  on  which  of 
these  dissociative  recombination  rate  coefficients  is  utilized  in  conjunction 
with  the  equilibrium  constant®^ 

Kg  =  (1.62  Tg  +  1.61  T^  +  2.2  T^)  x  10'^  exp(-2.8/Tg) 

Several  proces  -es  mcnl  ioned  in  the  report  have  no  measured  cross  sections, 

for  example,  the  disc:oc1i‘tion  of  0^  by  electron  impact.  In  such  cases,  the  colli- 
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Sion  rates  have  been  calculated  from  estimated  cross  sections,  which  may  in 
themselves  be  off  by  an  order  of  magnitude.  Other  rates  are  obtained  iy  analogy 
with  similar  processes  and  may  not  be  accurate  to  better  than  a  factor  of  2.  How¬ 
ever,  a  good  number  of  important  processes  have  been  measured  over  a  wide  tempera¬ 
ture  range,  and  are  reasonably  accurate  which  makes  the  description  of  ionizati(.i 
and  deionization  in  air  not  a  hopeless  task.  Indeed  we  expect  that  this  air 
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chemistry  code  can  predict  electron  densititos  in  reasonable  agreement  with  labor¬ 
atory  measurements  (a  factor  of  2  or  better).  Finally,  many  of  the  electron  im¬ 
pact  ionization  and  excitation  rate  coefficients  are  given  here  for  a  limited 
electron  temperature  range  for  convenience  only.  Their  values  at  higher  tempera¬ 
tures  are  also  available  and  could  be  obtained  from  us  upon  request. 

Laser  Preakdown  as  Predicted  bv  CHWIR 

The  calculated  response  of  the  ambient  atmospheric  air  (300K)  to  a  pulse 

Q 

of  COg-laser  radiation  is  shown  in  Figs.  k-S.  The  laser  intensity  is  5x10 

Watt/cm^  and  is  applied  for  x  =  0.13  usec.  It  should  be  noted  that  breakdown 

P 

9  2 

does  not  occur  for  laser  intensities  below  3x10  Wett/cn  ,  in  accord  with  exper- 

62 

imental  observations  ;  this  breakdown  threshold  is  tne  intensity  at  which  the 
loss  of  electrons  due  to  attachment,  exactly  balances  the  creation  of  electrons 
due  to  impact  ionization. 

Figure  4  plots  the  predicted  electron  density  The  two  curves  show  the 
different  results  obtained  using  two  different  dissociative  recombination  rates 

X  Q 

for  the  process  NO  +  e  ♦  N  +  0,  as  measured  by  Walls  and  Dunn  and  Huang, 

Biondi,  and  Johnson^®.  Results  presented  in  Figs.  5  and  6  are  based  on  the  dis¬ 
sociative  recombination  rates  measured  by  Biondi  and  his  colleagues. 

Figure  5  shows  the  evolution  of  the  positive  ions  N^^,  0*^  and  NO^. 

At  late  times,  the  dominant  ion  is  NO*,  which  balance  the  electron  density 

Figure  6  plots  N,  0,  and  NO.  The  reaction  N  +  0  ♦  NO*  +  e  and  its  inverse 
are  responsible  for  the  ionization  observed  at  late  times. 

Finally,  this  code  is  being  utilized  for  diagnostic  purposes  of  several  air 
ionization  experiments  at  NRL,  and  the  results  will  be  forthcoming  as  NRL  reports. 
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